Beyond-mean-field theories witii zero-range effective interactions. A way to handle 

the ultraviolet divergence 
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Zero-range effective interactions are commonly used in nuclear physics and in other domains to 
describe many-body systems within the mean-field model. If they are used within a beyond-mean- 
field fi'amework, contributions to the total energy that display an ultraviolet divergence are found. 
We propose a general strategy to regularize this divergence and we illustrate it in the case of the 
second-order corrections to the equation of state (EOS) of uniform symmetric matter. By setting 
a momentum cutoff A, we show that for every (physically meaningful) value of A it is possible to 
determine a new interaction such that the EOS with the second-order corrections reproduces the 
empirical EOS, with a fit of the same quality as that obtained at the mean-field level. 

PACS numbers: 21.60.Jz, 21.30.-x,21.65.Mn,02.60.Ed 



In many-body systems contact interactions are reason- 
able approximations to the realistic finite-range forces 
that can be employed in cases where the interaction range 
is smaller than the typical length scale represented by 
the interparticle distance. The main advantage of using 
zero-range interactions is that the equations to handle are 
greatly simplified. Two examples of commonly used con- 
tact interactions are the Skyrme forces which are quite 
popular in nuclear physics [1], and the contact interac- 
tions with coupling strengths depending on the s-wave 
scattering length which are employed for dilute atomic 
gases (see, e.g., Ref. [I])- In general, effective interactions 
contain parameters that must be adjusted to reproduce 
a given set of observables (for instance, in the nuclear 
case, binding energies and radii of a few selected nuclei, 
and bulk properties of nuclear matter). Usually this is 
done at the self-consistent mean-field level. When going 
beyond this level, it is likely that the effective interac- 
tion has to be redetermined. However, trying to discuss 
this issue is impossible in the case of contact interactions, 
since going beyond mean field, e.g., by including second- 
order corrections, implies dealing with contributions to 
the total energy which contain an ultraviolet divergence. 

Many authors have studied the perturbation series for 
the energy of the electron gas [3|-l5| and of nuclear sys- 
tems [6|-|9|. In all cases, the interaction is finite-range 
and no ultraviolet divergence appears. In this paper, we 
apply to a simple case a strategy to handle the contact in- 
teractions up to second order. We include a momentum 
cutoff A among the parameters of the interaction, and 
we show that for every value of A the other parameters 
can be determined in such a way that the total energy of 
the system with second-order contributions remains the 
same. This strategy, and the formulas we will show, are 
quite general and can be applied to different Fermi sys- 
tems. For the numerical application we restrict ourselves 
to symmetric nuclear matter, treated with a simplified 
zero-range interaction. 



The second-order terms that contribute to the total en- 
ergy in uniform matter and diverge in the case of contact 
forces are shown in the lower part of Fig. [TJ The diver- 
gence is caused by the integration on q and is somehow 
unphysical since the high-momentum states are certainly 
outside the scale at which effective forces are to be used. 

In the case of effective interactions between point- 
nucleons, the cutoff A must certainly be smaller than the 
momentum associated with the nucleon size, i.e., smaller 
than « 2 fm""'^. In fact, these interactions are used to de- 
scribe giant resonances or rotational bands of nuclei and 
consequently the scale should be even smaller, perhaps 
around 0.5 fm~^. However, our procedure is tailored on 
the basic idea of renormalizable quantum-field theories: 
so, in principle, for any value of A a new set of param- 
eters is found which leads to the same equation of state 
(EOS) -including the second-order contributions-. We 
show below that this is mathematically possible owing to 
the fact that the second-order correction is well-behaved 
as a function of the density. Thus, we imitate the QED 
idea that the bare electron mass and charge can be chosen 
for every different energy cutoff in such a way that the 
physical values of mass and charge are always obtained. 

An ultraviolet divergence appears already at the mean- 
field level in the Hartree-Fock-Bogoliubov (HFB) [13] 
or Bogoliubov-de Gennes (BdG) ,11] models when zero- 
range forces are employed in the pairing channel. In this 
case, sophisticated regularization schemes exist in which 
the irregular term of the pairing field is suppressed and 
the dependence on the energy cutoff is eliminated. These 
techniques are commonly adopted for example in atomic 
physics in BdG models |2|. |l2|. and have been also em- 
ployed in nuclear physics [l3| , but they cannot be directly 
applied to the case of interest studied in the present work. 

Let us write the zero-range force as 

V{v^,V2)^g5{v^-V2). (1) 
To make contact with the Skyrme interactions fH the 
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strength g is written as to + ^tap" and this corresponds 
to the so-called (toj^s) model that is a simplification 
of the usual Skyrme model where the spin-dependent, 
velocity-dependent and spin-orbit terms are dropped. If 
the Skyrme force is viewed as a G matrix (thus including 
ladder diagrams), the introduction of second-order con- 
tributions would in principle imply a double counting. 
However, our attitude is to consider as a matter of fact 
our interaction as phenomenological, and our framework 
as an effective theory where parameters are readjusted 
according to which diagrams are explicitly introduced. 

Normalizing the single-particle wave functions within a 
box of volume f2, the Hartree-Fock (HF) potential-energy 
contribution (upper part of Fig. [T]) is equal to 



Figure 1: First- and second-order diagrams for the total en- 
ergy in uniform matter. Labels refer to momentum states. 



This second-order correction is given by 
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where d = (n^ — n)/2, n being the level degeneracy (4 
in the case of symmetric nuclear matter), and v = ^. 
The energy per particle or EOS is obtained by adding 
the kinetic contribution and reads in our case 
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Eq. ^ coincides with the EOS obtained with the SkP 
parameter set , in which no contribution coming from 
the velocity-dependent terms appears and the effective 
mass coincides with the bare mass. Since SkP, as all 
the Skyrme sets, has been fitted to reproduce within HF 
the basic features of the nuclear EOS, one can consider 
its associated energy per particle as a benchmark which 
must be reproduced with reasonable accuracy for every 
value of A when the second-order correction is included. 



where C — — 6mi7^(27r) ^ft ^ for symmetric nuclear 
matter, and G{q) reads 



G{q) 



d^kid^k2 



(5) 



/fei,fe2<fcF,|ki-(-q|,k2-q|>fcF g2 -I- q • (ki - k2 
Some details about the evaluation of G(g) and AE with 
a zero-range force are recalled in the Appendix; finally, 
we write AE{p)/A as x{p) ^ /(p, oo), with 
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/(p, 00) = — uduFi{u) + J uduF2{u) ] (7) 
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where u — q/2kF. When the cutoff A is introduced, 
the last integral has A/2kp as upper limit and the corre- 
sponding quantity is denoted by /(p. A). The expressions 
for Fi (u) and F2 (u) are given in the Appendix. The an- 
alytical expression of /(p. A) is 
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In Fig. 2(a), E/A+AE/A is plotted for different values (solid black line). The correction AE/A is also shown in 
of A and compared with the SkP mean- field values of E/A Fig. 2(b). One observes that the second-order correction 
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Figure 2: (Color online) (a) E/A + IS.E/A as a function of 
the density and for different values of the cutoff A. The SkP 
mean-field EOS (solid black line) is shown for comparison, 
(b) Correction l^E/A for different values of A. 
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Figure 3: (a) Density-dependent x coefficient. Other panels: 
the values of /(p, A) for A = 0.5 (b), 1 (c) and f.5 (d) fm'V 



causes a shift of the saturation point to lower densities. 
The curves are calculated using the SkP parameters (cf. 
Table I). For A = 1.5 fm"^, the maximum correction is 
already comparable with the energy per particle at the 
saturation point, i.e., w 15 MeV. 

To better understand the behavior of /S.E/A, we plot 
x{p) and I{p,A) separately in Fig. |3l The ultraviolet 
divergence is visible at all the values of the density, by 
comparing the trends in panels (b), (c) and (d). At the 
same time, one observes also a divergent behavior of / 
for p that is dictated by the upper limit A/{2kp) 



of the second integral in Eq. ([7]). It can be seen in 
the analytical expression of Eq. ([8]). /(p. A) multiplies 
the coefficient x{p) which goes to zero like kp when the 
density goes to zero. One may worry that for large values 
of A, the divergent behavior of the integral dominates 
when /9 — i> 0, so that no regularization is possible for very 
low density. However, this happens for values of A which 
are physically meaningless, namely for A « 250-300 fm^^ . 

For each value of A we can perform a least square fit to 
determine a new parameter set SkPA, such that the EOS 
including the second-order correction matches rather well 
the one obtained with the original force SkP at the mean- 
field level. This is our main result and it is illustrated in 
Fig. 4(a). The refitted parameters are listed in Table I 
(together with the saturation point). In Fig. 4(b) we dis- 
play, for purely mathematical illustration, the refit done 
with the extreme value of A = 350 fm~^. 




0.04 0.08 0.12 0.16 0.2 0.24 0.28 0.32 
P(fm-') 



Figure 4: (Color online) (a) Second-order-corrected equations 
of state compared with the reference equation of state (SkP 
at mean-field level), (b) Extreme case of A = 350 fm~^. 



In summary we have shown that, for any physically 
meaningful value of the cutoff A it is possible to find a 
contact interaction that can be used in a mean-field plus 
second-order corrections framework and can describe sat- 
isfactorily the empirical EOS. The quality of the fits, that 
can be judged from Fig. 4 and the x^-values in Table I, 
demonstrates that the strategy we have outlined in a gen- 
eral fashion works, in practice, for the case of symmetric 
nuclear matter treated with a simplified contact force. 



One can foresee various applications to the studies of 
strongly correlated fermion systems in all the domains 
of many-body physics where zero-range forces are used 
and where beyond-mean-field theories are necessary for a 
more accurate treatment of complex correlations. In nu- 
clear physics, models beyond mean-field theories where 
ultraviolet divergences appear owing to the use of a zero- 



range interaction are: (i) models where multiparticle- 
multihole configurations are introduced [l5| . (ii) second 
random-phase-approximation models fl6| , and (iii) mod- 
els that take into account the coupling between single- 
particle degrees of freedom and collective vibrations |17l - 
l20l l. In these further development of our tech- 

nique may be envisaged. 
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Table I: From the second line, columns 2, 3 and 4: parameter sets obtained in the fits associated with different values of the 
cutoff A compared with the original set SkP (first line). In the fifth column the x^/-^- value (x^ divided by the number of fitted 
points) associated to each fit is shown. In columns 6 and 7 the saturation point is shown. 







to (MeV fm^) 


h (MeV fm^+^") 


a 


X'/N 


po (fm-«) 


E/A{po) (MeV) 




SkP 


-2931.70 


18708.97 


1/6 




0.16 


-15.95 


A = 


0.5 fm"^ 


-2352.900 


15379.861 


0.217 


0.00004 


0.16 


-15.96 


A = 


= 1 fm-i 


-1155.580 


9435.246 


0.572 


0.00142 


0.17 


-16.11 


A = 


1.5 fm"^ 


-754.131 


8278.251 


1.011 


0.00106 


0.17 


-16.09 


A = 


= 2 fm^^ 


-632.653 


5324.848 


0.886 


0.00192 


0.16 


-15.82 


A = 


350 fm"^ 


-64.904 


360.039 


0.425 


0.00042 


0.16 


-15.88 
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Appendix 

G{q) in Eq. ([5|) can be rewritten by expressing all wave 
vectors in units oikp'- 
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where the domains Di are D{ki) = {ki < 1, jk; + q| > 1}. 
By introducing y — aq, and the unit vector q = q/l^l, 
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1) First case: q > 2. In this case, |ki + q| > 1 and 
1^2 — q| > 1 are satisfied when ki < 1 and k2 < I. Eq. 
leads to 
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2) Second case: < g < 2. One can apply the tech- 
nique of Ref. 21], with this change of variables: 



dpfip,q) 
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With the change u = q/2 we obtain AE{p)/A as x{p) ^ 
/(p, 00) [cf. Eqs. dH) and with 
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where x = q - n and n is a unit vector n 
case, Eq. ((TU]) becomes 
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The energy correction can be written as 



AE = A-KCki 



q^dqv^Giiq) 



q^dqv^G2{q) (14) 
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